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ABSTRACT 

Although hclioscismology has been used as an effective tool for studying the physical mechanisms 
acting in most of the solar interior, the microscopic and dynamics of the deep core is still not well 
understood. Helioseismological anomalies may be partially resolved if the Sun captures light, non- 
annihilating dark matter particles, a currently discussed dark matter candidate that is motivated by 
recent direct detection limits. Once trapped, such particles (4-10 GeV) naturally fill the solar core. 
With the use of a well-defined stellar evolution code that takes into account an accurate description 
of the capture of dark matter particles by the Sun, we investigate the impact of such particles in its 
inner core. Even a relatively small amount of dark matter particles in the solar core will leave an 
imprint on the absolute frequency values of gravity modes, as well as the equidistant spacing between 
modes of the same degree. The period separation for gravity modes could reveal changes of up to 3% 
for annihilating dark matter and of up to 20% for non-annihilating dark matter. This effect is most 
pronounced in the case of the gravity dipole (I = 1) modes. 

Subject headings: elementary particles — dark matter — Sun: interior — Sun: helioseismology — 
stars: interiors 



1. INTRODUCTION 

Evidence for the existence of dark matter in the Uni- 
verse is well establis hed by cosmological observations 
(|Komatsu et al.l [2009f ). and its influence is also neces- 
sary in order to explain the f ormation of the struc ture of 
the current epoch Universe (jSpringel et al.l IM PS') . Var- 
ious studies suggest that dark matter is constituted of 
massive, neutral, weakly interacting non-baryonic parti- 
cles. Furthermore, independent considerations from par- 
ticle physics also pro pose the existence of such particles 
(jBertone et al.|[2005l ). 

Recently, the resul ts from several particle physics de- 
tection experiments (IBernabei et al.l 120 lCt lAhmed et all 
120101: lAalseth et all 120101: iFitzpatrick et all 120101 ) have 
been interpreted in terms of weakly interacting massive 
dark matter particles (WIMPs) with a mass smaller than 
10 GeV, and with spin-dependent elastic scat tering cross- 
sections with baryon s as la rge as 10 -32 cm 2 ([Angle et all 
2008; IBehnke et all 120081) 1 or spin- independent elas- 
tic scattering cro ss-section with ba ryons of the or- 
der of 10" 40 cm 2 (jAhmed et al.ll2009D . The relatively 
weak experimental limits on spin-dependent interac- 
tions of dark matter particles are especially of inter- 
est for the Sun with its large abundance of protons, 
for which both self-annihilating cold dark matter par- 
ticle scenarios (|Bergstrornl I2009D and non-annihilating 
(jFrandsen .V SarkaHl JO lTTT have previously been studied 



Investigation of the effects of the accumulation of dark 
matter particles by the Sun may therefore be an impor- 
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tant complement to direct detection searches for light 
WIMPs. These particles are trapped in the Sun's interior 
when they collide with nuclei and lose (linear) momen- 
tum, and drift into the Sun's core. Collisions of captured 
particles with the local baryons transfer and redistribute 
thermal energy, and lower the central temperature by 
a few percent. In thermal equilibrium, the kinetic en- 
ergy of dark matter p articles is balanced by the local 
gravitational potential (ISpergel fc PressH lQSSl An esti- 
mation of the radius of the dark matter core is given by 
r x ~ (9kT c / 4ttG p c m p ) yjrn p /m x where m p and m x are, 
respectively, the mass of the proton and the mass of the 
dark matter particle, T c and p c are the central tempera- 
ture and the central density of the Sun's core, and k and 
G are, respectively, the Boltzmann and Newton gravi- 
tational constants. This expression approximately gives 
the radius of the dark matter core in the Sun's interior. 
It follows that the more massive a dark matter particle, 
the smaller is the radius of the dark matter core, and 
the less important is the impact of dark matter in the 
evolution of the Sun. 

The Sun shows a complex pattern of surface oscilla- 
tions whose restoring forces are produced either by com- 
pressibility or buoyancy. The pressure perturbations give 
rise to acoustic sound waves in the high-frequency part 
of the spectrum, and buoyancy variations drive gravity 
waves in the low-frequency range of the spectrum. The 
small amplitude surface perturbations observed in the 
Sun can be described as a sum of eigenstates. Each 
eigenstate has a spatial counterpart that is defined by a 
spatial eigenfunction that depends on the thermodynam- 
ical structure of the background state (the Sun's internal 
structure), and a time-dependent eigenfunction that is 
characterized by the frequency v n ,i- The numbers I and 
n are positive integers, known as th e degree and radial 
order of the modefe.g.. lGouglj []"993'). 

During the last 50 years, accurate measurements of 
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frequency values were obtained for more than 5 thou- 
sand acoustic mode frequencies i^„. This achieve- 
ment is the result of the combined eff orts of several 
obser vational networ ks, such as BISON(|Chaplin et al.1 
[19961) and GONG (IHarvev et al.l fl996h . followed by 
the heliosei smic experiments from th e SoHO mis- 
sion: GOLF dGabriel et alJfl995l VIRGO dFrohlich et al.l 
[19951) and MDI-SOI (IScherrer et al.l Il995lh " Seismol- 
ogy has provided the most powerful tool to probe the 
Sun's interior to date. Accurate frequency measure- 
ments of large numbers acoustic modes have been made, 
including the radial and dipole global acoustic modes 
which penetrate more deeply into the core. Unfor- 
tunately, the diagnostics of the Sun's core provided 
by these acoustic modes are still insufficient. How- 
ever, the possible discove ry of gravity mod es by cur- 
rent ongoing experiments (|Garcia et al.l [2007) or future 
ones could be a breakthrough towards fully understand- 
ing the physics of the Sun's core and establishing the 
possible existence of dark matter. Th e impact of dark 
matter particles in the Sun's co re JCox et all |l990; 
IDearborn et~aTlll990t iKaplan et al.lll991[ ) and its conse- 
quences on the acous t ic spec trum were previously anal- 
ysed by lLopes et al.l (l2002bl) among other s, and more 
recent ly by iCumberbatch et al.l (|2010f ) and iTaoso et al.l 
(|201(fl . 

In this Letter, we have computed a series of solar mod- 
els evolving within dark matter halos for which we have 
estimated the impact of dark matter particles in the 
present-day Sun. As a reference, we use an updated solar 
standard model that provides a seismic diagnostic simi- 
lar to other solar standard models found in the literature. 
Furthermore, we have computed the gravity mode oscil- 
lation spectrum for such dark matter scenarios. 

2. HELIOSEISMOLOGY AND THE SUN'S EVOLUTION IN 
DARK MATTER HALOS 

The capture rate of dark matter particles in the 
Sun was computed in detail by taking into account 
their scattering wit h the differe n t chemical spe c ies of 
the solar plasma jGouldl Il987t iGondolo et "all 120041 : 
iCasanellas fc Lopesll2010l lbl). The presence of dark mat- 
ter inside the star can change its evolution by two mech- 
anis ms: (i) changing the transport of energy inside the 
star (ISpereel fc Pressl[l985l [Gould fc Raffeltlll990al) . or 
(ii) by contributing to the energy producti on that sus- 
tains the gravitati onal c ollapse of the star fllocco et al.l 
[200l iFreese et al.l [200l ICasanellas fc Lopesll2009D . In 
the case of the Sun , the latter mechanism is negligible 
(jLopes et al.ll2002bD . The efficiency of the energy trans- 
port by dark matter particles depends, among other fac- 
tors, on the scattering cross-sections of dark matter par- 
ticles on baryons. 

Dark matter particles might have a total scatter- 
ing cross-section a t with several components: spin- 
independent interactions with all chemical elements in 
the Sun, and spin-dependent interactions with hydro- 
gen. The presence of such particles contributes to 
the local energy transport. The efficiency of the heat 
transport by dark matter particles is dependent on the 
value at, relative to the fiducially critical cross-section 
a c = (m p /M)R 2 ~ 8 10 _36 cm 2 , where m p is the pro- 
ton mass and M and R are the mass and radius of 
the Sun. The optimal heat transport depends on the 
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Fig. 1. — Comparison of temperature profiles between the so- 
lar standard model and models of the Sun evolving in different 
dark matter halos. The halos are constituted of self-annihilating 
and non-annihilating massive particles with the following masses: 
(i) self-annihilating particles (coloured curves): 3 Gev (red), 4 
GeV (blue), 5 GeV (green), 6 GeV (cyan), 8 GeV (magenta), 10 
GeV (dashed-red) and 12 GeV (dashed-blue); (ii) non- annihilat- 
ing particles (black curves, from low to higher temperature): 5 
GeV, 7 GeV, 10 GeV and 50 GeV. The particles interact with 
baryons with a spin-dependent scattering cross-section of the or- 
der of 3 10 — 33 cm 2 (i: coloured curves) or 2 10 -35 cm 2 (ii: black 
curves), and spin-independent scattering cross-sections of the or- 
der of 10 — 40 cm 2 (negligible). The product of the self-annihilation 
cross-section and the relative velocity of colliding particles is of or- 
der of 3 10 _33 cm 3 s~ 1 (ixoloured curves) or 10 °cm 3 s — 1 (ii:black 
curves; suppressed annihilation). The black curve with the highest 
central temperature corresponds to the solar standard model. 

dark matter particle scattering cross-section cr t . For 
particles with small mean free paths (diffusion regime) 
or with large cross-sections (at > a c ), the conductiv- 
ity falls as a c /a t . Similarly, for particles with large 
mean free path (Knudsen regime) or lower cross-sections 
(a t < a c ), the conductivity falls as at/a c . Furthermore, 
a suppr ession factor is required in this latter transport 
regime ([Gould k Raffeltl[l990il lbh. Both regimes arc in- 
cluded in our stellar evolu tion code (|Lopes et al.l [2002b: 
ICasanellas fc Lop"esl 120091 ) . The numerical computation 
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of such effects on the transport of energy can be accom- 
plished by means of changing the production of energy 
or by changing the coefficient for radiative transfer. We 
have implemented both methods. The results obtained 
in both cases are similar. In this work, the transport 
of energy by dark matter is computed as a coefficient 
for the radiative transfer. This seems to be more re- 
alistic and is numerically more stable. Several models 
of evolution of the Sun were computed in different dark 
matter scenarios of suppressed and non-suppressed an- 
nihilation dark matter particles. Our reference model is 
an updated solar standard m odel (jTurck-Chieze fc Lopesl 
11993c lAsplund et al.l [2005) that shows an acoustic 
seismic diagnostic similar to other solar standard 
models (TTurck-Chieze et al.l l2004at iBahcall et all 120051: 
iGuzik fc Mussackll2010HTurck-Chieze et al.ll2oTo( ). We 
include self-annihilating and non-annihilating massive 
particles with the following masses: self-annihilating par- 
ticles from 3-12 GeV and non- annihilating particles 
from 5-50 GeV. The particles interact with baryons with 
a spin-dependent scattering cross-section 3 10~ 33 cm 2 
or 2 10~ 35 cm 2 , and a spin- independent scattering 
cross-sections of 10 _40 cm 2 (in fact, this is negligible). 
The product of the self-annihilation cross-section and 
the relative velocity of colliding particles taken to be 
3 10 _33 cm 3 s~ 1 or 10 _50 cm 3 s _1 (for the case of sup- 
pressed annihilation). 
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Fig. 2. — Comparison between the solar standard model and dif- 
ferent solar models evolved within self-annihilating dark matter 
halos: (a) sound speed radial profile and (b) density radial pro- 
file. The red dotted curve corresponds to the difference between 
inverted profiles and our solar standard model (see text). See cap- 
tion of Figure 1 for the details on self-annihilating scenarios (case 
i: coloured curves). 
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Fig. 3. — Comparison period table P; n between the solar stan- 
dard model and models of the Sun evolving in different dark mat- 
ter halos (scenario i). The periods P; „ shown corresponds to the 
dipole gravity modes (I = I). This corresponds to dark mat- 
ter halos constituted of annihilating massive particles, (case i: 
coloured curves). See caption of Figure 1 for the details of the 
self-annihilating scenario. 

In the case of non-annihilating particles, the numerical 
model shows a 2% reduction of the central temperature of 
the Sun's core due to the lightest massive particles (Fig- 
ure 1). This temperat ure reduction seems to be stronger 
than that obtained bv lTaoso et al.l ((2010) . Additionally, 
the square of the sound speed and density profiles incur a 
variation of between 0.8% and 2% (Figure 2). In the case 
of non-annihilating particles, the decrease in the central 
temperature is of the order of 8 % for a particle with a 
mass of 12 GeV, and even larger for lighter particles. 

Figure 2 shows the profiles of sound speed and den- 
sity in the case of annihilating dark matter scenar- 
ios compared with the inverted profiles. The inversion 
of the sound speed and density profiles were done us- 
ing th e seismic data of the BISON and GONG net- 
works (jBasu et~aT1[20C]9T) . This seismic data is consis- 
tent with the previous high accuracy measurements done 
by th e GOLF and MDI instru ments of the SoHO mis- 
sion (|Turck-Chieze et al.ll200ll) . It is evident that the 
present acoustic seismology is not able to probe the inner 
core of the Sun accurately. In summary, our knowledge 
about the Sun's interior is quite accurate up to 10% of 
the S un's radius (jTurck-Chieze et al.H2004bt IGarci'a et al.l 
2007). The dark matter is expected to accumulate in 
the inner core of the Sun. An accurate description of 
such deep layers of the Sun's interior can be obtained if 
gravity modes are observed. 

3. GRAVITY MODES AND THE EXISTENCE OF A DARK 
MATTER CORE 

An additional probe of the Sun's interior, capable of 
finding small dark matter effects occurring in the Sun's 
core, can be performed through gravity modes. We 
have computed the gravity modes of oscillation of the 
present Sun for different da rk matter halos scenarios 
(jChristensen-D alsgaard 2008J). The changes in the spec- 
trum of gravity modes, due to the existence of dark mat- 
ter in the Sun's core are visible even in the case of very 
small amounts of dark matter accumulated in its core. 
In the case of self-annihilating dark matter scenarios, 
the period for gravity dipole (I = 1) mode changes are 
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from a few tens of seconds up to 10 minutes (see Fig- 
ure 3). In non-annihilating dark matter scenarios, the 
effect is much stronger. For very light particles, the grav- 
ity dipole modes could have their periods reduced by 80 
minutes. This existence of a dark matter core also affects 
the period spacing of gravity modes, a quantity identical 
to the large separation for acoustic modes. In principle, 
this should provide the first clear indication of the ex- 
istence of a gravity oscillation pattern in the observed 
spectrum. A qualitative expression for the large period 
separation can be obtained for the case of gravity modes 
with low degree I and high-order n modes, where the 
period Pt, n (= uf*) is g iven by 



Pl,n = 



Po 



VUT+Tj 



I 



o 



1 



Pl,nJ ' 



with 



Po = 2tt 2 



\N\ 



dr 



(1) 



(2) 



where N is the buoyancy, and is a phase term sensi- 
tive to the lay ers below the base of the convective zone 
(]TassouH 11980) . This expression, and in particular P a , 
tells us that gravity mode frequencies are determined by 
the density stratification of the core, through the buoy- 
ancy. The value of P is of the order of 34.10 minutes 
in the case of our solar standard model. Usually, the 
value of P is computed from the large period separa- 
tion, §Pi tn = Pi t n — Pi, n -i' In the case of dipole modes 
(I = 1), it reads 5P\ tn — P Q /y/2. In the case when the 
dark matter core becomes isothermal, the value of P a 
is strongly affected. The structure differences between 
the solar standard model (ssm) and different dark mat- 
ter scenarios (dm) can be estimated by computing large 
separation period differences, which is equivalent to mea- 
sure period differences, AP^ n — Pf^—Pf^ 71 (I and n fix). 
It follows that APi^/P?™ 1 AP /P SS "\ The structure 
differences produced by the presence of dark matter in 
the Sun's core leads to a significant change in the period 
separation. Figure 4 shows the large period separation 
for several self-annihilating dark matter scenarios in the 
case of gravity dipole modes (I — 1). The period spacing 
in such models can experiment changes of up to 3 %. In 
the case of non-annihilating scenarios the period spacing 
can be reduced by as much as 20 %. 

4. CONCLUSION 

The possible existence of a dark matter core, even in 
the case of non-annihilating dark matter scenarios of low 
mass particles, will be always contained within the first 5- 
10% of the solar radius. Such a dark matter core will very 
likely have a small temperature gradient where dark mat- 
ter particles are present and, in the case of strongly inter- 
acting dark matter scenarios, it will form a fully isother- 
mal core. Such an isothermal core will produce two dis- 
tinct signatures in the gravity modes: (i) by changing the 
frequency values of gravity modes, (ii) and by changing 
the large period separation between gravity modes of the 
same degree and consecutive radial order. 

The impact of the dark matter in the evolution of the 
Sun and its impact in the helioseismology acoustic data 
as well as in the solar neutrinos observables has been 
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Fig. 4. — The large separation for gravity modes: The large sep- 
aration was computed from periods Pi, n , SPi n = P; n — P; n _i 

Pol + 1) shown in Figure 3. This was computed for dipole 
gravity modes (I = 1). These dark matter halos are constituted 
of self-annihilating massive particles (case i: coloured curves). See 
caption of Figure 1 for details of the self-annihilating scenario. 



addressed by severa l authors(e.g. , iDearborn et alJl!990b 
Kaplan et all 11991b lLopes et all l2002bt Hopes fe Silkl 



2002t Hopes et alj|2002al lLopes fc Si lk 2010). More re- 



cently a specific study has been done to explore a new 
class of fundamental particle candidates which are able to 
prod uce significant chang es in the structure of core of the 
Sun (iCumberbatch et al1l20Tot iFrandsen fe Sarkarll2010t 
iTaoso et al.1 l2010f T In this study, we have explored in 
more detail how the gravity spectrum could be modified 
by the presence of such types of dark matter particles in 
the Sun's core. 

Nevertheless, other physical processes participating in 
the evolution of stars need to be better understood 
in order to take full account of the effects caused by 
dark matter in the Sun's core. In particular, the new 
CNO composition has lead to a smaller central tem- 
perat ure than the one required by neutrino detectio n 
(e.g., iGuzik fc Mussacldl2010t iTurck-Chieze et alJlMoh . 
The inclusion of well-known physical processes in the 
Sun's evolution, such as differential rotation, meridional 
circulation, magnetic breaking, formation and evolu- 
tion of the solar tachocline layer and sjalar dynamo (e.g., 
iCharbonne au 2005; iPassos fc Lopesl l2008f ) , as well as 
the transport of angular momentum by gravity waves 
and/or magnetic fields (among other dynamical pro- 
cesses) will lead to minor structur e changes through- 
out t he evolution of the star (e.g., ITurck-Chieze et all 
2010). In some cases, such dynamical processes will in- 
crease the discrepancy between solar neutrino measure- 
ments and acoustic seismology. In other cases, hydrody- 
namical processes will lead to reductions of the (small) 
differ ences between the solar mod el and observations 
(e.g. JGaraud fc Bodenh cimcr 20l3). Nevertheless, a full 
quantitative account of the physical processes involved 
in the evolution of the star is fundamental to finding the 
signature of dark matter in the Sun's core. 

The possible discovery of gravity modes by current 
space missions such as SoHO or the new generation ex- 
periments SDO and PICARD will be of major impor- 
tance in the search for dark matter inside the Sun. In 
addition to these, the GOLF-NG instrument is specifi- 
cally designed for a future solar mission. The goal of this 
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instrument is to probe the very central region of the Sun, 
around 0.5% of the solar radius, the most likely place to 
find dark matter. 

We thank the anonymous referee for advice in im- 
proving this paper. We gratefully acknowledge the au- 
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